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Abstract 

Due to the inereasing number of industrial applieations of eleetro-aetive polymers 
(EAPs), there is a growing need for eleetromeehanieal models which accurately capture 
their behavior. To this end, we compare the predicted behavior of EAPs undergoing ho¬ 
mogenous deformations according to three electromechanical models. The first model is 
a continuum based model composed of the mechanical Gent model and a linear relation¬ 
ship between the electric field and the polarization. The electrical and the mechanical re¬ 
sponses according to the second model are based on the polymer microstructure, whereas 
the third model incorporates a neo-Hookean mechanical response and a microstructural 
based long-chains model for the electrical behavior. In the microstructural motivated mod¬ 
els the integration from the microscopic to the macroscopic levels is accomplished by the 
micro-sphere technique. Pour types of homogenous boundary conditions are considered 
and the behaviors determined according to the three models are compared. The differ¬ 
ences between the predictions of the models are discussed, highlighting the need for an 
in-depth investigation of the relations between the structure and the behaviors of the EAPs 
at microscopic level and their overall macroscopic response. 
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1 Introduction 

Dieleetrie elastomers (DEs) are materials that deform under electrostatic excitation. Due 
to their light weight, flexibility and availability these materials can be used in a wide variety 
of applications such as artificial muscles Q, energy-harvesting devices [|2l|3l, micropumps 
IIH, and tunable wave guides [I51[6l|, among others. The principle of the actuation is based on 
the attraction between two oppositely charged electrodes attached to the faces of a thin soft 
elastomer sheet. Due to Poisson’s effect, the sheet expands in the transverse direction. Toupin 
ITTI , in his theoretical work, found that this electromechanical coupling is characterized by a 
quadratic dependence on the applied electric field. This was later verified experimentally by 
Kofod et al. (HI. However, DEs have a low energy-density in comparison with other actuators 
such as piezoelectrics and shape memory alloys (fl. Eurthermore, their feasibility is limited 
due to the high electric fields (~ 100MV/m) required for a meaningful actuation as a result of 
the relatively low ratio between the dielectric and elastic moduli (91 HOl . Specifically, common 
flexible polymers have low dielectric moduli while polymers with high dielectric moduli are 
usually stiff. Nevertheless, a few recent works suggest that this ratio may be improved. Huang 
et al. (TTlI demonstrated experimentally that organic composite EAPs (electro-active polymers) 
experience more than 8% actuation strain in response to an activation field of 20 MV/m. The 
experimental work of Stoyanov et al. lfT2ll showed that the actuation can be dramatically im¬ 
proved by embedding conducting particles in a soft polymer. In parallel, theoretical works 
dealing with the enhancement of coupling in composites also hint at the possibility of improved 
actuation with an appropriate adjustment of their microstructure (51 11^ [T4l [15111^ . 

The above findings motivate an in-depth multiscale analysis of the electromechanical cou¬ 
pling in elastic dielectrics which is inherent from their microstructure. In this work we consider 
the class of polymer dielectrics. A polymer is a hierarchical structure of polymer chains each 
of which is a long string of repeating monomers. We start by analyzing the behavior of a single 
monomer in a chain. Next, the response of a chain is obtained by a first level integration from 
the single-monomer level to the chain level. Einally, the macroscopic behavior of the polymer 
is obtained by a higher level summation over all chains. In this work, we utilize existing con¬ 
stitutive models for the chains and concentrate on the higher level summation from the chain 
level to the macroscopic-continuum level. To this end the physically motivated micro-sphere 
technique, that enables to extend one dimensional models to three dimensional models by ap¬ 
propriate integration over the orientation space, is exploited (TTIIT^ . Accordingly, this method 
lends itself to the characterization of polymer networks since the single chain is often treated 
as a 1-D object which is aligned along the chain’s end-to-end vector lfT9ll20ll . 

The response of a polymer subjected to purely mechanical loadings was extensively inves¬ 
tigated at all length scales. A macroscopic level analysis and models describing the behavior 
of soft materials undergoing large deformation, such as polymers, were developed by Ogden 




Noy Cohen, Andreas Menzel and Gal deBotton 


3 


m. The microscopic level analysis of Kuhn and Griin [l22l yielded a Langevin based eon- 
stitutive relation and paved the way to various multiseale models sueh as the 3-ehain model 
(231, the tetrahedral model (24l|23| and the 8-ehain model (261 . Corresponding micro-sphere 
implementation of the Langevin model was carried out by Miehe et al. [|T9l . 

The electric response of dielectrics to electrostatic excitation was examined by Tiersten (27l 
and Hutter et al. (28], among others, macroscopically as well as through their microstructure. 
Starting with the examination of a single charge under an electric field, the relations between 
the different electric macroscopic quantities, such as the electric displacement, the electric field 
and the polarization, and microscopic quantities, such as the free and bound charge densities 
and the dipoles were defined and analyzed. 

The study of the response of dielectrics to a coupled electromechanical loading initiated 
with the pioneering work of Toupin (Til, who performed a theoretical analysis at the macro¬ 
scopic level. Later on, an invariant-based representation for the constitutive behavior of EAPs 
was introduced by Dorfmann and Ogden (291 . Subsequently, Ask et al. (3011311 . Gei et al. (51 
and Jimenez and McMeeking (3^ investigated the possible influence of the deformation and 
its rate on the electromechanical coupling. Thylander et al. (20l made use of a corresponding 
micro-sphere technique at the chain level. By employing macroscopic constitutive models for 
the mechanical and the electrical behavior of the polymer chains, a few boundary value prob¬ 
lems were solved by means of the numerical implementation of the micro-sphere technique and 
the finite element method. Initial multiscale analyses of the electromechanical response were 
performed by Cohen and deBotton (3311341 . The present work focuses on the implementation 
of different electromechanical models to EAPs experiencing homogenous deformations under 
various types of boundary conditions and examination of their predicted response. 

We begin this work with the detailed description of the different macroscopic and micro¬ 
scopic models and the presentation of the micro-sphere framework. Next, the micro-sphere 
technique is used to compute the macroscopic behavior of dielectrics with randomly oriented 
and uniformly distributed dipoles experiencing the macroscopic rotation. In section 4 we de¬ 
termine and compare the behavior of a polymer according to three electromechanical models 
under different boundary conditions. The conclusions are gathered in section 5. 

2 Theoretical background 

Consider the deformation of a hyperelastic dielectric body subjected to electro-mechanical 
loading from a referential configuration to a current one. In the reference configuration, the 
body occupies a region Vq G with a boundary dVo. The referential location of a material 
point is X. In the current configuration, the body occupies a region E G with a boundary 
dV, and we denote the location of a material point by x. The mapping of positions of material 
points from the reference to the current configurations is x = ^(X), and the corresponding 
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deformation gradient is 

F = Vx9, (1) 

where the operation Vx denotes the gradient with respeet to the referential eoordinate system. 
The right and left Cauchy-Green strain tensors are C = F^F and b = FF^. The ratio between 
the volumes of an infinitesimal element in the current and the reference configurations, J = 
detF, is strictly positive. In the case of incompressible materials, which are of interest in the 
present work, we have 7=1. 

The induced electric field E on the body satisfies the governing equation 


VxXE = 0, 


( 2 ) 


in the entire space, where Vx denotes the gradient with respect to the current coordinate system. 
Consequently, we can define a scalar field, the electric potential <^, such that E = —Vx^- The 
electric displacement field is 

D = CqE 4" P, (3) 

where Cq is the permittivity of the vacuum and P is the polarization, or the electric dipole- 
density. We recall that in vacuum P = 0. In the absence of free charges the electric displacement 
field is governed by the local equation 


Vx-D = 0. (4) 

In the work of Il29l . the referential counterparts E^*^) and of the electric field and the 
electric displacement were determined. Specifically, 


e(0) = 

(5) 

D(o) =7F ^D, 

(6) 


where Vx x E*^®^ = 0 and Vx = 0. We note that unlike E*^®) and the referential 
polarization is not uniquely defined. In order to ensure that the referential polarization is energy 
conjugate to the referential electric field such that ■ p(®^ = E • P, we adapt the definition 

p(0)=7F^ip. (7) 

In accordance with our assumption that the dielectric solid can be treated as a hyper-elastic 
material, its constitutive behavior can be characterized in terms of a scalar electrical enthalpy 
per unit volume function W. We further assume that W can be decomposed into a mechanical 
and a coupled contributions, i.e. W (F, E) = Wq (F) 4- Wc (F, E), where Wq (F) characterizes the 
material response in the absence of electric excitation and Wc (F, E) accounts for the difference 
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between W with and without eleetrie exeitation [1^ l34l [35l |36l . Aeeordingly, the polarization 
is determined via 


P 


J dE ’ 


( 8 ) 


and the stress developing in the material can be written as the sum 


G — 7C-\-Gy., 


(9) 


where 


G/n 


1 dW (F) 

7 <9F 



is the mechanical stress due to the deformation of the material, 


( 10 ) 


;r = E(8)P, 


( 11 ) 


is the polarization stress stemming from the applied electric field in the dielectric, and 


Gy 


£o 


E E 


[E-ElI 


( 12 ) 


is the Maxwell stress in vacuum, where I is the second order identity tensor [l35l . We empha¬ 
size that this decomposition is purely modelling-based as in an experiment the total stress can 
be measured, but the contributions of the individual components cannot be distinguished. In 
this work, we consider incompressible materials that undergo homogenous deformations and 
therefore a pressure like term pi, which is determined from the boundary conditions, is added 
to the total stress. Assuming no body forces, the stress satisfies the local equilibrium equation 


Vx a = 0. 


(13) 


The electrical boundary conditions are given in terms of either the electric potential or the 
charge per unit area pa, such that D ■ n = —pa, where n is the outward pointing unit normal. 
Practically, in EAPs, pa is the charge on the electrodes. The mechanical boundary conditions 
are given in terms of the displacement or the mechanical traction t. Due to the presence of 
the electric field the stress in the vacuum outside of the body does not vanish. Therefore, 
the mechanical traction at the boundary is [a — C7v.] -6 = 1, where the expression for Gy, the 
Maxwell stress outside the material, is given in Eq. ( [T^ in terms of the electric field in the 
vacuum. 

2.1 Existing models for the behavior of dielectrics 

Within the framework of finite deformation elasticity the simplest constitutive model is the 
well-known neo-Hookean model that requires only one material parameter. The corresponding 
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strain energy-density funetion (SEDF) is 

Wo^^iF) = ^[h-3], (14) 

where /i is the shear modulus and /i = tr(C) is the first invariant of the right Cauehy-Green 
strain tensor. This model does not capture the lock-up effect observed in experiments and 
corresponds to a significant stiffening of the material at large strains. Gent OTl proposed a 
phenomenological constitutive model in which this effect is accounted for. The SEDF for this 
model 

= -^In (^1 - , (15) 

depends on two parameters, p and Jm- The latter is the lock-up parameter such that -|- 3 is 
the value of Ii at the lock-up stretch. Thus, the expression in Eq. ( [T5] ) becomes unbounded at 
h =Jm + 3 which captures this phenomenon. 

With regard to the response of the dielectric to electrical excitation a quadratic dependence 
of Wc on E is commonly assumed, leading to the linear relation 

P = ZE, (16) 

where x is the susceptibility of the material OSll . Consequently, the electric displacement is 

D = eE, (17) 

where e = Co -f is the permittivity of the material. We note that this linear relation is in 
agreement with the invariant based representation of Dorfmann and Ogden [[29l . Furthermore, 
an experiment carried out by Di Eillo et al. ll^ on VHB 4910 showed that this assumption is 
fairly accurate. 

Recent experiments with various types of polymers imply that the permittivity, and there¬ 
fore the relation between the polarization and the electric field, is deformation dependent 
ll40l l4Tl l42l l43l . A possible explanation for this dependency of the susceptibility on the defor¬ 
mation is related to the alteration of the inner structure of the polymer [|33ll34ll . 

In the polymer, the monomers in the chain can move or rotate relative to their neighbors thus 
providing the chains with a freedom to deform [l44l . In order to better understand the response 
of polymers to an electro-mechanical loading, their microstructure should be accounted for. 
This can be accomplished in terms of a multiscale analysis consisting of three stages: the first 
involves the examination of the behavior of the monomers, the second includes analysis of the 
response of the chains, and the third deals with the polymer behavior at the continuum level. 

Treloar [|45ll and Flory and Rehner [|24ll carried out a multiscale analysis of a polymer sub¬ 
jected to mechanical loading. It was assumed that the directions of the monomers, or links. 
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eomposing a ehain are random, and eonsequently it was found that the ehains are distributed 
aeeording to a Gaussian distribution. Based on statistieal eonsiderations and the laws of ther- 
modynamies, the variation in the entropy of the chain due to its deformation was determined. 
The overall variation in the entropy of the polymer is computed by summing the entropies of 
the chains. Remarkably, their result recovered the macroscopic neo-Hookean behavior. Fur¬ 
thermore, a comparison between the micro and macro analyses related the macroscopic shear 
modulus to the number of chains per unit volume Nq. Specifically, it was found that 

^i=kTNo, (18) 


where k and T are the Boltzmann constant and the absolute temperature, respectively. 

Due to the assumptions that lead to the use of Gaussian statistics, the lock-up effect was not 
captured in the above mentioned analysis. A more rigorous examination of the polymer behav¬ 
ior by Kuhn and Griin [1221 revealed that this phenomenon is a result of the finite extensibility 
of the chains. According to this analysis the SEDF associated with a polymer chain is 


W^^ = kTni 


( — 1 - -bln 




nilJ nil 


sinhf^-if^ 


(19) 


where I is the length of a link, n/ is the number of links in a chain, r is the distance between the 
two ends of the chain, and (•) is the inverse of the Langevin function 


i?(/5)^coth(/5)-f = 

P nil 


( 20 ) 


Assuming that all chains undergo the macroscopic deformation, i.e. r = Fro where r and ro are 
the current and referential end-to-end vectors, respectively, the stress associated with a chain is 


derived from Eq. (191, 


Gjf = kT^i'^^-^ ( — ) Fro®Fro, 
r \nil' 


( 21 ) 


where ro is a unit vector in the direction of ro and Eq. ([T0|) is used. Here, 


rQ = I y/nu 


( 22 ) 


is the average length of the referential end-to-end vectors [|44l 14611 . The quantity de- 

scribes the ratio between the end-to-end and the contour lengths of the chain, and from Eq. 


(221 it follows that 


— = VFf 0 ■ Fro ^ = VFf 0 ■ Ff o 

nil nil 


1 


(23) 
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It ean be shown that the limit 1 results in capturing the exper¬ 

imentally observed lock-up phenomenon. The lock-up stretch is associated with the chain 
undergoing the largest extension such that the stretch ratio of its end-to-end vector is 


^max — x/^l- 


(24) 


It is important to note that the first term in the Taylor series expansion of Eq. (21) about ^ = 0 
reproduces the Gaussian model [|22l |46l . 

A few works proposed models that consider specific finite networks of chains. The 3-chain 
model by Wang and Guth [123]| examines a network of 3 chains which are located along the axis 
of the principal directions of the deformation gradient. Flory and Rehner [|24l and Treloar Il251l 
proposed a network of four chains that are linked together at the center of a regular tetrahedron, 
and their other ends are located at the vertices of the tetrahedron. The tetrahedron deforms ac¬ 
cording to the macroscopic deformation while the chains experience different stretches. In the 
model proposed by Arruda and Boyce 8 representative chains in specific directions rela¬ 
tive to the principal system of the macroscopic deformation gradient are used to determine the 
macroscopic behavior. An anisotropic worm-like chain model in which no inherent alignment 
between the chosen and the principal coordinate systems is assumed was considered by Kuhl 
et al. [|47l . 

A multiscale level analysis of the response of polymers with Gaussian distribution to electro¬ 
mechanical loading was carried out by Cohen and deBotton [133]|. In this study, the changes in 
the magnitudes of the dipolar monomers due to the applied electric field and their rearrange¬ 
ment due to the mechanical deformation were accounted for. Following a model described in 
Stockmayer [|48ll . Cohen and deBotton 041 considered the class of uniaxial dipoles in which 
the dipole is aligned with the line segment between the two contact points of a monomer to its 
neighbors. Thus, taking | to be the unit vector along this line segment the dipole moment of a 
monomer is 


HIm = 




E, 


(25) 


where .fi is a material constant. The dipole of a chain composed of uniaxial dipoles with an 
end-to-end vector in the direction r is 


m. 




+ 

1_ 

r 

[I-3r®r] 


[ndl\ 



E. 


If we assume that all chains undergo the macroscopic deformation, then r = 
accordance with a second type of dipoles discussed in Stockmayer 


Frp 


(26) 


In 


V [Frol-[Frol 

, Cohen and deBotton 

proposed an expression for a transversely isotropic (TI) dipole, where the dipole is aligned 
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with the projeetion of the eleetrie field on the plane perpendieular to ^, 


m. 


2 




E. 


(27) 


An expression for the dipole of a ehain made out of transversely isotropie dipoles is determined 
by following the steps followed in the derivation of Eq. ( [2^ for the uniaxial dipoles. 

The resulting polarization of the polymer is determined by summing the dipoles of the 
ehain in a representative volume element of a volume via [|38]| 

(28) 

and the resulting polarization stress is eomputed via Eq. 

2.2 The micro-sphere technique 


Consider a unit sphere whose surfaee represents the directions of the referential end-to-end 
vectors. The directional averaging of a quantity • over the unit sphere can be approximated 
with a discrete summation 




• dA 


I 

/=1 




(29) 


where the index / = 1,..., m refers to a unit direction vector Tq where • ' is the value of quantity 
• in the direction and is an appropriate non-negative weight function constrained by 


= 1 [fT7l[T9ll^ . In general, Eq. ( [29| ) can be combined with a more general anisotropic 
distribution function [ISOllSTI . 

In the case of polymers the vectors ^ represent the directions of the end-to-end vectors of 
the polymer chains or, from a numerical point of view, the integration directions in orientation 
space. Eor randomly and isotropically oriented chains, or rather isotropic integration schemes, 
the vectors satisfy 


dO, 


Tq'W 


= 0 , 


/=1 


and 






r^w 


1=1 




(30) 


(31) 


In view of Eqs. (301 and ( pT) the micro-sphere technique naturally lends itself to the calculation 
of the macroscopic polarization and stress. Specifically, for a polymer with chains composed 
of dipoles and Nq chains per unit referential volume, Eq. may be written as 
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and the maeroseopie stress aeeording to Eq. (21) as 


m J \ / ' 


(33) 


where we use the notation suggested in Eq. ( [29| ). 

Bazant and Oh ifTTl demonstrated that a speeifie ehoiee of 42 direetions guarantees suffi- 
eient aeeuraey for the applieation diseussed in their work. We follow this eonjeeture where the 
integration directions and the corresponding weight functions are given in Table 1 of Bazant 
and Oh ifTTIl . We note that other integration schemes are available, as demonstrated by Waffen- 
schmidt et al. [[52]| . Ostwald et al. [|5^ and the references cited therein. 


3 Dielectrics with randomly distributed monomers 


Consider a model of a dielectric composed of hq monomers per unit referential volume, 
which are treated as mechanical rods and electric dipoles. The dielectric is subjected to a 
mechanical deformation, locally represented by F, and an electric field E. We assume that the 
electric field induced on a monomer by its neighbors is small in comparison with the applied 
electric field [[34ll . We examine first a dielectric with uniaxial monomers, the behavior of which 


is governed by the quadratic form in Eq. (251. If all of the dipoles experience the macroscopic 


rotation, i.e. ^=R^o where R = FC is a proper rotational tensor, then the polarization 
according to Eq. ([32]) is 


42 


P = no (me) (8)|o^ = ^^E, 

(=1 


UqA 


(34) 


where Eqs. (29) and ([3T|) are used. The corresponding polarization stress is 


;r = E(8)P= '^EGjE. 


(35) 


In the case of a dielectric composed of hq TI dipolar monomers per unit referential volume, 
cf. Eq. ( [27] ), which mechanically act as rigid rods, the same assumptions that led to Eq. 0 
lead to 


^ 42 

P = no (me) = «o 2 L 


1=1 






w^E = 




E. 


(36) 


). We note that the 
are identical to the 

exact expressions obtained by Cohen and deBotton ll34l . 


Accordingly, the expression for the polarization stress is given in Eq. (35 
polarization and polarization stress calculated in Eqs. (|34|), (|3^ and ([35|) 
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4 Dielectric elastomers 

We examine the behaviors of incompressible dielectric elastomers according to three dif¬ 
ferent models under various homogenous electromechanical loading conditions and compare 
between their predicted responses. To facilitate the comparison we assume that in the limit 
of infinitesimal deformations and small electric excitations all three models admit the same 
behavior. Specifically we assume that the initial shear modulus and electric susceptibility are 
/r = 0.1 MPa and % = SEq. In those models in which the lock-up effect is accounted for we 
choose the model parameters such that under purely mechanical biaxial loading the lock-up 
stretch is = 5. The precise models and the numerical values assumed for their parameters 
are as follows: 


1. The macroscopic model - the mechanical behavior is characterized by the Gent model 
( [T5] ) with the aforementioned shear modulus and = 47. The electric behavior is deter¬ 
mined according to the linear model ([T7]) with the initial permittivity £ = 4£q. 


The microscopic model - the Langevin model pT] ) is utilized in order to describe the me¬ 
chanical behavior, where = 25 is chosen to fit the assumed lock-up stretch according 


to Eq. p4|) and Nq = We employ the long-chains model p6l) with chains that are 


composed of = 100 uniaxial dipoles (Eq. 25) to characterize the dielectric response 
of the polymer. The material constant ^ is chosen such that = 3 Cq to ensure that 
the referential polarization is identical to the one admitted by the macroscopic model. 


3. The Gaussian model - the neo-Hookean model (14) with (18) are used, where Nq = 
in conjunction with the long-chains model ( |2^ to characterize the mechanical and the 
electrical behaviors, respectively. We assume that a chain is composed of = 100 
uniaxial dipoles (Eq. [^, where the chosen long-chains model constant is identical to 
the one determined for the microscopic model. 


In the following, we examine a thin layer of a polymer whose opposite faces are covered with 
flexible electrodes with negligible stiffness. The electrodes are charged with opposite charges 
so that the difference in the electric potential induces an electric field across the layer. Erom a 
mechanical point of view we consider four boundary conditions. In the first two cases differ¬ 
ent displacements are prescribed at the boundary and consequently the deformation gradient is 
defined. In the following two representative cases we set the traction on the boundaries. We 
choose a cartesian coordinate system in which the referential electric field is aligned with the y- 
axis and calculate the macroscopic polarization according to the microscopic and the Gaussian 
models via Eq. ( |^ . The polarization stress is computed according to Eq. G3- The me¬ 
chanical stress is computed via Eqs. ( [T5| ), ( [T4] ) and ( |2T] ) for the macroscopic, the Gaussian and 
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the mieroseopie models, respeetively. The pressure term is determined from the traetion free 
boundaries to whieh the eleetrodes are attaehed and subsequently the total stress is eomputed 
via Eq. Q. 

For eonvenienee we define the dimensionless normal stress along the x-axis d = ■ CTx 

and the dimensionless referential eleetrie field and referential eleetrie displaeement along the 
y-direction ■ y and ■ y, respeetively. In the following examples 

the current configuration counterparts of and are E = E-yandD=^D-y, 
respectively, as follows from Eqs. Q and 


4.1 Equibiaxial stretching perpendicular to the electric field 

In this case the material is stretched along the axes x and z such that Ax = = A. As 

stated previously, the y-axis is aligned with the referential electric field and due to the assumed 
incompressibility Xy = This setting is common in various experiments with EAPs 


Figs. §>) and ([^) depict o and D as functions of A^ and E, respectively. The curves 
with the squared marks correspond to the macroscopic model, the curves with the hollow cir¬ 
cle marks to the Gaussian model, and the curves with the filled circle marks to the microscopic 
model. The applied referential electric field is = 50 We point out that at A^ = 1 the di¬ 
mensionless stress according to the three models is not zero but very small. Fig. ([^) illustrates 
the stress increase at the lock-up stretch according to the macroscopic and the microscopic 
models. As expected, this effect is not observed when the Gaussian model is employed. Fur¬ 
thermore, since the electric field tends to stretch the material in the transverse plane, as long 
as the prescribed stretch is smaller than the electrically induced stretch, the overall stress is 
compressive. We emphasize that since the potential is held fixed, as the layer is stretched the 
current electric field increases and hence also the electromechanic ally induced stress. This 
gives rise to different types of loss of stability phenomena which are outside the scope of the 
current work. The reader is referred to the works by, e.g., Rudykh et al. llTSll . Dorfmann and 
Ogden Il54l . Bertoldi and Gei [|55ll and Shmuel et al. [0. Only when the prescribed stretches 
are large enough, the total stress becomes tensile. 

In Fig. ([^) we observe a linear dependence of the electric displacement on the electric 
field according to the macroscopic model, as follows from Eq. (17) and the assumed constant 
permittivity. Since the graph is plotted in terms of the dimensionless quantities its slope is unity. 
In contrast, the Gaussian and the microscopic models predict a stronger than linear increase in 
the electric displacement as we stretch the material. This is a result of the predicted increase 
in the permittivity due to the stretching of a polymer with chains made up of uniaxial dipoles 
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Figure 1: The dimensionless stress d (a) and eleetrie displaeement D (b) versus the streteh of 
the transverse plane and the dimensionless eleetrie field aeeording to the maeroseopie model 
(the eurve with the square marks), the Gaussian model (the eurve with the hollow eirele marks) 
and the mieroseopie model (the eurve with the filled eirele marks). 


4.2 Pure shear deformation in the plane of the electric field 

Once again we analyze a thin layer of a polymer whose opposite faces are covered with 
flexible electrodes with negligible stiffness. As before the y-axis is aligned with the electric 
field, but in this case the deformation of the material along the z-axis is constraint such that 
= 1. The layer is stretched along the x-axis such that Xx = X and the incompressibility 
condition yields Xy = j. 

Fig. ([^) depicts the dimensionless normal stress that develops along the x-axis as a func¬ 
tion of the stretch according to the macroscopic model (the curve with the square marks), 
the Gaussian model (the curve with the hollow circle marks) and the microscopic model (the 
curve with the filled circle marks) under the applied referential electric field = 100^. 
The inability of the neo-Hookean model to capture the lock-up stretch is again clearly de¬ 
picted. We also notice that there is a difference between the lock-up stretches predicted by 
the macroscopic and the microscopic models. This is because the lock-up stretch according to 
the Langevin model is determined by the maximum eigenvalue of the deformation gradient as 
seen from Eq. ( [2^ , whereas according to the Gent model it depends on the first invariant of the 
right Cauchy-Green strain tensor. Treloar [|46l presented experimental results demonstrating 
that polymers lock-up at different values under different types of deformations, and therefore 
we conclude that in this aspect the Gent model may be a better predictor. We wish to point out 
that the lock-up stretch according to the microscopically motivated 8-chain model of Arruda 
and Boyce [|^ . in which the chain behaves according to Eq. ( [^ , depends on /i as well and 
is able to capture the different lock-up stretch values under various states of deformation. A 
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Figure 2: The dimensionless stress d (a) and eleetric displaeement D (b) versus the axial 
streteh and the dimensionless eleetrie field aceording to the maeroseopie model (the eurve 
with the square marks), the Gaussian model (the eurve with the hollow eirele marks) and the 
microseopie model (the curve with the filled circle marks). 



comparison of different models and their calibration according to the data reported by Treloar 
||4^ was carried out by Steinmann et al. [f5^ . 

The curves with the square, hollow and filled circle marks in Fig. (|^) correspond to the 
macroscopic model, the Gaussian model and the microscopic model, respectively. Here, the 
dependency of D on £ is illustrated. Due to the constant permittivity, we again note the linear 
dependency predicted by the macroscopic model. The Gaussian and the microscopic models, 
which are based on the electric long-chains model, predict a change in the permittivity as a 
result of the mechanical stretch, in agreement with the experimental findings of Choi et al. 
[|40ll . Wissler and Mazza iHTTl . McKay et al. [l42l . Qiang et al. PGl . Since the deformation is 
dictated by the boundary condition, the Gaussian and the microscopic models predict the same 
electric behavior. 

4.3 Equibiaxial actuation normal to the electric field 


We once again examine a thin layer of a polymer whose opposite faces are covered with 
flexible electrodes with negligible stiffness. This time, however, the circumferential boundary 
of the layer is traction free, thus allowing the layer to expand in the plane transverse to the di¬ 
rection of the electric field in response to the electric excitation. We choose the same system of 


axes as in subsection 4.1 and, thanks to the symmetry of the loading, the deformation gradient 
is diagonal with = A. Due to the assumed incompressibility we have Ay = ^. 

Fig. ([^) displays the dimensionless referential electric field as a function of the in¬ 
duced stretch of the transverse plane according to the macroscopic model (the curve with the 
square marks), the Gaussian model (the curve with the hollow circle marks) and the micro¬ 
scopic model (the curve with the filled circle marks). The loss of stability discussed in the 
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Figure 3: The dimensionless referential electrie field (a) and eleetric displaeement 
(b) versus the stretch of the transverse plane and the dimensionless referential electric field 
according to the macroscopic model (the curve with the square marks), the Gaussian model 
(the curve with the hollow circle marks) and the microscopic model (the curve with the filled 
circle marks). 


works of Dorfmann and Ogden ll54ll . Bertoldi and Gei [[55l . Rudykh and deBotton ll57l and 
Shmuel et al. m is demonstrated again. We note that after the peak at ^ 0.7, even though 
the current electric field increases monotonically, the Gaussian model predicts a decrease in 
the referential electric field with an increase of the stretch. In an experiment where the referen¬ 
tial electric field is controlled, the macroscopic and microscopic models predict a jump in the 
planar stretch. This effect of a transition between two states is known as snap-through (I5^ |4|). 

The curves with the square, hollow and filled circle marks in Fig. (|^) correspond to the 
macroscopic, the Gaussian and the microscopic models, where the predicted dependencies of 
oiO) 

on in the direction of the electric field are depicted. Essentially, this plot illustrates 
the amount of charge per unit referential surface area as a function of the potential difference 
divided by the initial thickness of the layer. Initially, we observe an increase of the surface 
charge with an increase of the electric potential. However, beyond the peak at E^^"> ^ 0.7 
there is a reversed trend where, at equilibrium, the electric potential drops while the surface 
charge increases. This occurs in conjunction with the uncontrollable increase in the area of 
the actuator as shown in Fig. ([^). From a practical viewpoint this implies that beyond the 
peak, in a manner reminiscent of an electrical short-circuit, excessive current flows from the 
system electric source while the electric potential drops. We stress that due to the thinning of 
the layer the current electric field increases and may result in a failure of the DE due to electric 
breakdown. We also note that even though the same electric model is used in both the Gaussian 
and the microscopic models, there is a difference in the relations between the electric field and 
the electric displacement. This is due to the different mechanical deformations resulting from 
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Figure 4: The dimensionless referential eleetrie field (a) and eleetrie displacement (b) 
versus the stretch and the dimensionless referential electric field according to the macroscopic 
model (the curve with the square marks), the Gaussian model (the curve with the hollow circle 
marks) and the microscopic model (the curve with the filled circle marks). 


the applied electric field according to the two different models. 

4.4 Uniaxial actuation normal to the electric field 


We consider a setting reminiscent of the one considered in subsection 4.3, but this time 
the layer is free to expand only along the x-direction. Consequently, the deformation gradient 
components are Aj. = 1, = A, and Xy=j. 

Fig. (|^) shows the dimensionless referential electric field as a function of the stretch 
according to the macroscopic model (the curve with the square marks), the Gaussian model 
(the curve with the hollow circle marks) and the microscopic model (the curve with the filled 
circle marks). As mentioned previously, the predicted lock up stretches are A^“ = 5 and A^“ ~ 7 
according to the microscopic and the macroscopic models, respectively. Unlike the biaxial 


case described in subsection 4.3 this loading does not admit loss of stability according to the 
macroscopic and the microscopic models. The Gaussian model, that is based on the mechanical 
neo-Hookean model, reaches a peak at 1 and then monotonically decreases. 

Fig. (|^) depicts the dependence of on According to the Gaussian model, repre¬ 
sented by the curve with the hollow circle marks, no significant changes in the surface charge 
are observed as we initially increase the electric potential difference. However, beyond the peak 
of 1, this model predicts an unstable behavior as a result of the electrical long-chains 

model. In order to maintain equilibrium according to the macroscopic and the microscopic 
models an increase in the charge on the electrodes requires an increase in the voltage differ¬ 
ence between them. Thus, the Gaussian model is admitting a behavior that is qualitatively 
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different from the behaviors of the other two models. 

5 Concluding remarks 

We determined the behavior of an incompressible polymer undergoing homogenous defor¬ 
mations according to three electromechanical models under four types of boundary conditions. 
The first model incorporates well-known macroscopically motivated constitutive relations for 
the mechanical and the electrical behaviors. The second, microscopic model, combines me¬ 
chanical and electrical models stemming from the microstructure of the polymer. The third 
model assumes a Gaussian distribution of the polymer chains and accordingly the mechani¬ 
cal and the electrical behaviors are determined. We comment that the material parameters nj 
and ni, which denote the number of dipoles and links, respectively, are used as fitting parame¬ 
ters and therefore the electrical long-chains model and the Langevin model are not consistent. 
Further investigations in this regard is needed. 

In the first two representative examples we apply a referential electric field by setting the 
potential difference between the electrodes and controlling the deformation. In the following 
two examples we apply a referential electric field and set the traction on the boundaries of the 
polymer. In order to determine the polarization and the stress resulting according to the micro¬ 
scopic models we make use of the micro-sphere technique. A comparison between the results 
shows that the macroscopically and the microscopically motivated models predict different be¬ 
haviors. Therefore, this work encourages a further and a more rigorous investigation into the 
connection between the two scales aimed at deepening our understanding of the micro-macro 
relations and the mechanisms which control the actuation. Moreover, analysis of this type may 
open the path to the design and manufacturing of polymers with microstructures that enable to 
improve the electromechanical coupling. 
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